Solution of electric-field-driven tight-binding lattice in contact with fermion reservoirs 
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Electrons in tight-binding lattice driven by DC electric field dissipate their energy through on-site 
fermionic thermostats. Due to the translational invariance in the transport direction, the problem 
can be block-diagonalized. We solve this time-dependent quadratic problem and demonstrate that 
the problem has an oscillatory steady-state. The steady-state occupation number shows that the 
Fermi surface disappears for any damping from the thermostats and any finite electric field. Despite 
the lack of momentum scattering, the conductivity takes the same form as the semi-classical Ohmic 
expression from the relaxation-time approximation. 

PACS numbers: 71.27.+a, 71.10.Fd, 71.45. Gm 



I. INTRODUCTION 

Nonequilibrium phenomena in lattice is one of the old- 
est and most fundamental problems in solid state physics. 
In conventional solids, the acceleration due to external 
field is relatively small compared to electronic energy, 
and various scattering mechanisms make the transport 
diffusive enough so that the small field approximation 
has often been used. The quantum Boltzmann method 
has been applied effectively^! and linear response limit 
has been widely used in the literature. However, recent 
progress in nano-devices and optical lattice systems has 
made rigorous high-field formalism necessary to under- 
stand their non-perturbative effects such as the Bloch os- 
cillation. In such regime, understanding the interplay of 
non-perturbative field-effect and the many-body physics 
has emerged as one of the most pressing problems in 
nano-science. 

Combining the nonequilibrium and quantum many- 
body effects is an extremely challenging task. Much ef- 
fort has been exerted towards understanding strong cor- 
relation physics in quantum dot physics, especially the 
nonequilibrium Kondo problem. Analytical theories 3 ^ 
and many numerical methods have been proposed along 
the time-dependeniP^, anc j steady-state simulation d 10 | 11 l 
In such systems with localized interacting region, the im- 
portant question of energy dissipation could have been 
side-stepped, and the existence of steady-state has not 
been a major issue. 

In the past few years, non-perturbative inclusion of 
electric-field and many-body effects has been one of the 
important issues in the field. Theor ies for lattice nonequi- 
librium have been formulate d 12 * 13 !, mostly based on the 
dynamical mean-field theory (DMFT) for an s-orbital 
tight-binding (TB) lattice with on-site interaction^^. 
Although a long-held belief in solid-state transport has 
been that, under a finite electric-field, the Fermi sea is 
perturbatively shifted by a drift velocity, many calcula- 
tions performed under the DMFT framework have sug- 
gested that the system approaches a steady-state with 
infinitely hot electron gas even for small field. With in- 
clusion of proper dissipation mechanism, correct models 
should expect the Boltzmann picture of displaced Fermi 



surface at small fields and a recovery of the Bloch oscil- 
lation in the high-field limit. 

However, it is unclear so far which of the approx- 
imations, such as single-band approximation without 
Landau-Zener tunneling or the nature of on-site inter- 
action, are responsible for the rather peculiar long-time 
states. One of the goals of this paper is that we provide 
exact solutions to one of the dissipation models with on- 
site fermion thermostats and give analytic understanding 
of the problem, and guide possible future modeling. 

Due to the nature of the one-body reservoirs, the prob- 
lem can be solved exactly (see Fig. [IJ. With identical 
reservoirs on each site, the Hamiltonian can be block- 
diagonalized according to the wave- vector of electrons in 
the transport direction. The block-diagonal Hamiltonian 
can then b e exac tly solved by a time-dependent perturba- 
tion theorjff^EiD U gi n g the nonequilibrium Green function 
theory. The calculation of the occupation number sup- 
ports the numerical observations that the Fermi surface 
of the unperturbed equilibrium is completely lost with 
effectively infinite temperature, in contrast to the con- 
ventional picture of the displaced Fermi sea under finite 
electric field. Despite the lack of momentum scattering, 
the DC electric current of this model recovers the familiar 
semi-classical Boltzmann equation resultP^. 

II. MODEL 

We study a quadratic model of a one-dimensional s- 
orbital tight-binding model connected to fermionic reser- 
voirs (see Fig. [IJ under a uniform electric field E. The 
effect of the electric field is absorbed in the temporal 
gauge as the Peierls phase ip(t) — eEat to the hopping in- 
tegral The time-dependent Hamiltonian then reads 
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with d\ as the (spinless) electron operator on the tight- 
binding chain on site i, c\ a with the reservoir fermion 
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FIG. 1: One-dimensional tight-binding lattice of orbital di 
under an electric field E. Each lattice site is connected to an 
identical fermionic bath of {da} with the continuum index a 
along the reservoir chain direction. 



states connected to the site i with the continuum index 
a along each reservoir chain. Here we do not specify 
the explicit connectivity of the reservoir chains, but each 
chains are assumed to have an identical dispersion re- 
lation e a . Notice that the electric field is applied only 
on the tight- binding chain {d\}. The coupling between 
the TB site and the reservoir is given by the identical 
tunneling parameter g. The Peierls phase tp(t) is given 
as 



¥>(*) = 



0, for t < 
fit, for t > 



(2) 



fi = eEa is the Bloch-oscillation frequency due to the 
electric field. 

We note that the whole system has discrete trans- 
lational symmetry in the transport direction and the 
Hamiltonian is readily block-diagonalized with respect 



to the wave- vector k as d\ = y N d 1 e lkRj dj and 

C L = \f^Y,j &ikRic \ a ( with latticc sitcs R j = aj 
and the number of sites along the TB chain Nd), 



Hit) = E 



-2 7 cos(fc + ip{t))d\d k + E 



(3) 



Here e<j(fe) = — 2~fcos(k) is the tight-binding dispersion 
at zero E-field. Then each /c-sector can be treated and 
solved separately. So from now on, we suppress the fc- 
subscript until necessary with the following Hamiltonian, 

H k (t) = -2~fcos(k + (p(t))dU + J2e a cic a 

(X 

-gJ2( c U + H.c). (4) 



It is important to note that the fc-dependence enters the 
problem as k + <p(t) for t > 0. This problem is simply a 
resonant level modePSl where the level is modulated si- 
nusoidally for t > 0. Considering these two observations, 
we can infer that different k only adds different phase to 
k + (p(t) and steady-state properties such as the average 
occupation number Ukif) = (d* (t)c(t)) are expected to be 
k- independent. This shows that the on-site fermion ther- 
mostat problem completely obliterates the Fermi surface 
at any non-zero field and, even in the small damping 
limit, and does not restore the conventional picture of 
shifted Fermi surface due to an external field. 



III. SOLUTION FOR OCCUPATION NUMBER 
AND CURRENT 

The time-dependent Hamiltonian Q can be exactly 
solved by the nonequilibrium Keldysh Green function 
method^]. We write the Hamiltonian as Hk{t) — 
Hq + V(t) with the time-independent unperturbed part 
Hq = .Hfc(O) and the time-dependent perturbation as 
V{t) = H k (t) - H k (0), 

V(t) = -2 7 [cos(fc + tp(t)) - cos(fc)] dU = v(t)dU. (5) 

When the perturbation is one-body on discrete states the 
lesser and greater part of the self-energy is zero, and the 
lesser d-Green function G K is expressed only in terms of 
the transient term, symbolically written as^ 



[I + G r V]G^ [I + VG a ] 



(6) 



and the retarded Green function G r is given by the usual 
Dyson's equation 



G r = G r + G r VG r , 



(7) 



where the matrix product means convolution-integrals in 
time. 

First with the retarded functions, the non-interacting 
limit has the time-translational symmetry and 



G r (t-t') = -i6(t-t') 



de 2 , r2 e 



e(t-t') 



= - 2 0(i-t')e" l£d(fe)(t ^*' ) ~ r| *~ t ' 1 , (8) 

where we we use a flat-band DOS for the reservoir in the 
infinite-band limit with the hybridization broadening T = 
Trg 2 N(0) and the density of states of the fermion bath 
^(0) = EJitod- Writing G r (t,t') = G r (t - t')g r (t,t'), 
Eq. ^ becomes 



g r (t,t') = l-i J dsv(s)g r ( S ,t'), 



which can be solved as 

g r (t,t') = exp 



i J v(s)ds 



(9) 



(10) 
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FIG. 2: Expectation value of occupation number rik(t) of 
wave-vector k and the current Jk(t) for 7 = 1, Q = 0.5, 
T = 0.1 and at k = tt/2 + 0.1. After the initial time of r -1 , 
transient behavior dies out and the expectation values reach 
steady oscillation. 



and finally we have for the full retarded Green function 
G r {t,t') = -^(t-t')e"' ied(fc)(t - t ' ) ~ r| *~ t ' l exp -if v(s)ds 

For the same-time argument for G < we have the 
Dyson's equation 

G<{t.t) = G$(t,t) 

+ f [G r (t,s)w(s)G<(s,<) + G<(t,s)w(s)G Q (s,i)]ds 





+ 11 G r (t,s)v(s)G^(s,s / )v{s / )G a (s' ,t)dsds' '.(12) 
Jo Jo 

We set the initial lesser Green function with the half-filled 
reservoir at zero temperature as 



G<(t,i') = i [ duj- ^. 2 , r2 e- 



j(t-t') 



(13) 



After some straightforward steps, the occupation number 
for the wave-vector k, n k (t) — —iG < (t,t), becomes 



da) 



T/tt 



(Lo-e d {k)f+T^ 



(14) 



1-i / dsv(s)e l( "- £d(fc)+ir)(t - s) -^> (s ' )<ls ' 



n k (t) = 



Fig. [2] shows the above n k (t) numerically evaluated for 
7 = 1, = 0.5, T = 0.1 and at k = tt/2 + 0.1. Due to 
the exponential factor e~ r ('~ s \ the integral converges to 
a steady-state oscillation state after time t ~ and 
the transient behavior dies out. Therefore, for long-time 
behavior, the time-integral range [0, t] can be changed to 



[— 00, t] for easier analytic treatment. After an integral- 
by-parts and some straightforward steps, we have 



n k (t) = 



duj x 



^ g e ~i(LU+ir)s-i(2-//(l) sin(fc+fi(t+s)) 



(15) 



An identity for Bessel functions J n (x) 

00 



Ax cos 6 



(16) 
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can be used to perform the integrals as 



n k (t) 



E 



T (2j_\ j ( 2j\ i(m-n)(k+Ot) 
<Jn\ q )<Jm\ q ;c 

-(m-n)Q + 2iT 

2 n 2 + r 2 



1 TO 

2 ° g n 2 ^ 2 + T 2 



with 



Xr, 



7r + tan 



tan 



n57 



(17) 



(18) 



As pointed out earlier, it is important to note that the 
k and t-dependence is only in the exponential factor as 
e i(m-n)(fc+sit)_ Therefore, taking a time-average over the 
period 2iri}~ 1 gives /c-independent results, which have 
been confirmed in numerical calculations. For instance, 
the DC limit of n k (t) only involves the terms of to = n 
in Eq. (17 1, and we have 



1 OO 

5 £ 



-i 2 



(19) 



For electronic models without thermostats, it has been 
argued that the effective electron temperature quickly 
becomes infinity^ from ^-independent occupation num- 
ber fife. The same conclusion applies to this model with 
on-site fermionic thermostats. 

Now we turn to the calculation of electric current, 

Mt) = nt) n k (t) = 2 7 sm(fe + Qt)n k (t). (20) 

Due to the sine-function, the DC current has contribu- 
tions only from to — n = ±1 in Eq. (17 1. After some 
manipulations, we have 



Jk — 



2 7 r 



7r(o 2 + 4r 2 ) 



riog 



(to - i) 2 ^ 2 + r 2 



J i, 



(21) 



As in the case for n k (t), the DC limit of J k (t) becomes 
independent of k. The total current is shown in Figs. [3] 
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FIG. 4: Contour plot of DC current as a function of damping 
and electric field. 



FIG. 3: DC current as a function of the damping V and elec- 
tric field Q = eEa. For small field, the current has a linear 
dependence on the i5-field showing an Ohm's law-like behav- 
ior. As E increases, the Bloch oscillation behavior takes over 
and the DC current decreases. The dashed lines are the sim- 



plified expression, Eq. ( 22 1 



and [4] Similar plot has been obtained in the interact- 
ing model from numerical calculation of Hubbard model 
connected to fermion batifSD 

It is instructive to simplify the above expression in the 
limit of f2, r <C 7 where the DC current is reduced to the 
expression 



J : 



4 7 m 



7r(fF + 4r 2 )' 



(22) 



Detailed derivation is provided in Appendix. This ap- 
proximate expression is shown as dashed lines in Fig. [3] 
Despite that the formula was derived for !!,r « 7, it 
shows remarkable accuracy to the DC current for a wide 
range of V and E. 

It is also interesting to note that a similar formula 
has been derived for a super-lattice system with Ohmic 
scattering within the semi-classical Boltzmann transport 
equatiorPn Although the current has the same depen- 
dence on the damping and the electric field, it should be 
emphasized that the two models have quite different scat- 
tering mechanism where in the Boltzmann approach^! 
the momentum relaxation is explicitly built-in while in 
our case the lattice wave-vector scattering does not hap- 
pen and a very different Fermi surface structure results. 
In the low-field limit, the current (22 1 recovers the form 



of the Drude conductivity per electron, 



7 r2 Et 

-p oc — , 
7rl m 



with 7 ~ 1/m* and T ~ 1/r. 



(23) 



IV. CONCLUSIONS 

Calculations on electron transport with fermionic ther- 
mostat confirm salient features of numerical results, 
such as uniform occupation number for all fc-vectors in 
the Brillouin zone and the Ohmic-like limit of electric 
current. In particular, the simplified electric current, 
Eq. (22), recovered the semi-classical Boltzmann trans- 



port result. However, the nature of the oscillatory steady- 
state is quite different from the conventional solid-state 
transport picture. In the presence of a static electric field, 
this model predicts that the Fermi surface completely dis- 
appears for any finite damping T and electric field E after 
the time-evolution of the interval At > max(f2 , T -1 ). 
The applicability of this model may be justified for stud- 
ies on short-time dynamics and dissipation. However, us- 
ing this model to study the strongly correlated nonequi- 
librium steady-state requires much caution since the uni- 
form on-site coupling to reservoirs does not seem to pre- 
serve the Fermi surface. 

We may speculate what could be the minimal elements 
for a realistic model which may bridge the conventional 
picture of Fermi surface shift and the non-perturbative 
field-effect, hence eventually accessing the nonequilib- 
rium strongly correlated steady-state. One may consider 
momentum scattering (i) in phononic Ohmic bath (ii) in 
fermionic bath with site-disorder and (iii) Landau-Zener 
tunneling in multi-band dispersion. 
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Appendix A: Derivation of current at small field Using the asymptotic expression^ for x/D,,j/Q — > oo, 



For both r, f2 <C 7, we expand Eq. ( 21 ) to the leading 
order of m as, 



2 7 r 



ir(n 2 + 4r 2 ) 

2mr^ 2 



Jm—\ 

_ 1 mfl 



202 , P2 ^ 2 ^ tan r 



(Al) 



Rearranging the summation and using the identity 
x[J m _i(x) + J m+ i(x)} = 2mJ m (x), we write 



■h 



2 7 r 



tt(^ 2 +4F 2 ) 

mm 2 



[Jm-1 + ^m+l] X 



, h f2 tan 

m 2 Q 2 + r 2 r 



L °° J l0 (W>2 7 ) 



the integral simplifies to 



2-) 



-27 vr ^47 2 - x 2 V^ 2 + r 2 



x tan — ax. 

r, 



(A2) 



2m 



n{n 2 + 4r 2 ) 



/ 777,rr2 mf2\ ^ n nm ^ r<7, the second term in the parenthesis 

EmfJJ™ ( — — —7 + tan -1 — — ) . dominates and we have 
m \m 2 tt 2 + T 2 T J 



For f2 7, we define x = m£l in the regime m = x/fl 3> 
1, the summation becomes 



xr _•■ x\ dx 



2- 



xdx 



4jm 



n(n 2 +4r 2 ) J v /4 7 2 - a; 2 7r(fi 2 + 4r 2 )- 



(A3) 
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